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Multifractal analysis of the metal-insulator transition in the 3D Anderson model I: 

Symmetry relation under typical averaging 
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The multifractality of the critical eigenstate at the metal to insulator transition (MIT) in the three- 
dimensional Anderson model of localization is characterized by its associated singularity spectrum 
f{a). Recent works in ID and 2D critical systems have suggested an exact symmetry relation 
in /(a). Here we show the validity of the symmetry at the Anderson MIT with high numerical 
accuracy and for very large system sizes. We discuss the necessary statistical analysis that supports 
this conclusion. We have obtained the f{a) from the box- and system-size scaling of the typical 
average of the generalized inverse participation ratios. We show that the best symmetry in f{a) 
for typical averaging is achieved by system-size scaling, following a strategy that emphasizes using 
larger system sizes even if this necessitates fewer disorder realizations. 
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I. INTRODUCTION 

The Anderson model of localization has been a sub- 
ject of intense analyticali*^!^ and numerical^i^ studies 
for decades. Anderson in his seminal paper— had demon- 
strated that, at absolute zero temperature and in the 
absence of external fields and electron-electron interac- 
tions, a sufficiently strong disorder can drive a transition 
from a metallic to an insulating state (MIT). The scaling 
theory of localization^ has shown that such a transition 
arises generically for systems with dimension d > 2 . 

A characteristic feature of this critical transi- 
tion is the strong multifractality of its wavefunction 
amplitudes i^i^i^ d°i^^ The critical eigenstate being neither 
extended nor localized reveals large fluctuations of wave- 
function amplitudes at all length scales. The character- 
ization of multifractality is most often given in terms of 
the singularity spectrum /(a). It can be computed from 
the g-th moments of the inverse participation ratio Pq 
which defines the scaling behaviour of Pq ex A'^*^'^"^-'^"'^' 
with length XJ^ The anomalous exponents Aq determine 
the scale dependence of the wave function correlationsi^ 
and separate the critical point from the metallic phase 
(for which Aq = 0). By carrying out such a multifrac- 
tal analysisi^iiiii^ (MFA) various critical properties can 
be obtained such as the critical disorder Wc,— the posi- 
tion of the mobility edges and the disorder-energy phase 
diagram,— as well as the critical exponents of the local- 
ization length^iiSiSa 

From an analytical viewpoint not much is known about 
the singularity spectrum. An approximate expression can 
be obtained in the regime of weak multifractality, i.e. 
when the critical point is close to a metallic behaviour. 
This applies to the Anderson transition in d = 2 -I- e di- 
mensions with e <C 1. In this case a parabolic dependence 
of the spectrum is found as f{a) ~ d~[a ~ {d + e)] /4e)^ 
which in turn implies Aq 2± —eq{q — 1). Although the 
parabolic approximation has turned out to be exact for 
some models^ its validity, in particular for the integer 



quantum Hall transition, is currently under an intense 
debate^^i^ due to the implications that this result has 
upon the critical theories describing the transition. 

It is only in the thermodynamic limit where a true 
critical point exists and hence the true critical /(a) can 
be obtained. Since the numerical characterization of the 
multifractal properties of \ip\'^ at the MIT can only be 
obtained from finite-size states, one therefore has to con- 
sider averages over different realizations of the disorder. 
Due to the nature of the distribution of Pq ^^^ one would 
normally take the typical average which is exactly the ge- 
ometric average of the moments of Pq over all contribu- 
tions. The use of typical averaging for the MFA has been 
successfully implemented in various studiesi^ d^d^d"^ 

Remarkably, it was recently argued that an exact sym,- 
metry relation should hold for the anomalous scaling 
exponents 
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A, = Ai_„ (1) 

which for the singularity spectrum is translated into. 



f{2d-a) ^ f{a)+d-a. 



(2) 



This relation implies that the singularity strength a must 
be contained in the interval [0, 2d] and that the values of 
/(a) for a < d can be mapped to the values for a > d, 
and vice versa. We note that the parabolic /(a)Si is 
in perfect agreement with this form of the singularity 
spectrum provided that /(a) is indeed terminated at 
and 2d. Numerical calculations have since then sup- 
ported this symmetry in /(a) in the one-dimensional 
power-law random-banded-matrix model^l and the two- 
dimensional Anderson transition in the spin-orbit sym- 
metry class j^^i^^ In the present work wc numerically ver- 
ify that this symmetry in the singularity spectrum also 
holds in the three-dimensional (3D) Anderson model. In 
order to address this hypothesis with sufficient accuracy, 
we have considered the box- and system-size scaling of 
the typical average of Pq in computing the /(a). We 
discuss which numerical strategy will produce the best 



possible agreement with the symmetry and we highhght 
the statistical analysis that must be used to observe the 
reported symmetries with sufficient confidence. In a re- 
lated publication, we also address this problem using the 
ensem6/e-averaged approachS^ and the reader may wish 
to compare both articles for a more complete picture of 
the MFA at the Anderson transition. 



II. THE MODEL AND ITS NUMERICAL 
DIAGONALIZATION 

We use the tight-binding Anderson Hamiltonian in lat- 
tice site basis as given by 



n = Y,e, \i){i\ + Y,U, |i)0-|, 



(3) 
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where site i = (a;, y, z) is the position of an electron in a 
cubic lattice of volume V = L^ , Uj are nearest- neighbour 
hopping amplitudes and £,; is the random site potential 
energy. We consider £,; to have a box probability distri- 
bution in the interval [— Wc/2, Wc/2], where Wc is taken 
to be the strength of the critical value of the disorder. 
We assume Wc = 16.5, above which all eigenstates arc 
localised.i^Si^iHi^ Furthermore the hopping amplitude is 
taken to be i = 1 and periodic boundary conditions are 
used to minimize boundary effects. 

The L^ X L^ Hamiltonian is diagonalized using the 
JADAMILU package^! which is a Jacobi-Davidson im- 
plementation with an integrated solver based on the 
incomplcte-iC/-factorization package ILUPACKj ^'^i^^ We 
have considered eigenstates tp = X^iV'iK) o^ly i^^ the 
vicinity of the band centre E = Q, taking about five eigen- 
states in a small energy window aX, E — Q for any given 
realization of disorder. A list of the number of states and 
the size of -0^ used for each L is given in Table HI For 
computing the singularity spectrum using the so-called 
box-size scaling approach (see section IIV|) . the largest 
system size we used is L = 240 with 95 eigenstates. This 
translates into 1.3 x 10^ values of wave function ampli- 
tudes ipi. For the system-size scaling, we used all system 
sizes in TableHl A critical eigenstate for L = 240 is shown 
in Fig. m 



III. MULTIFRACTAL ANALYSIS 

A. Basic definitions 

Let I'f/'ip be the value at the z-th site of a normalized 
electronic wavefunction in a discretized d-dimensional 
system with volume U^. If we cover the system with Ni 
boxes of linear size I, the probability to find the electron 
in the /c-th box is simply given by 




FIG. 1; Multifractal eigenstate for the 3D Anderson model 
at _E = and Wc = 16.5 for linear system size Z/=240 with 
periodic boundary conditions. The 410075 sites with proba- 
bility IV'jP twice larger than the average l/L^ are shown as 
bo xes w ith volume \ipj\^L^. The 26097 boxes with IV'il^i^ > 
2\/1000 are plotted with black edges. The color/grey scale 
distinguishes between different slices of the system along the 
axis into the page. 



The /Xfc(Z) constitutes a normalized measure for which we 
can define the q-ih moment as 



N, 



ps) = Y.i'i^^)- 



(5) 



k=l 



The moments Pq can be considered as the generalized 
inverse-participation ratios (gIPR) for the integrated 
measure /ifc(/), reducing to the wave function itself in 
the case / = 1 (in units of the lattice spacing) and to the 
usual IPR P2 for q = 2. The general assumption underly- 
ing multifractality is that within a certain range of values 
for the ratio A = l/L, the moments Pq show a power-law 
behaviour indicating the absence of length scales in the 
system^ 
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The mass exponent T{q) is defined as 



r{q) 
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for metals. 
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at the MIT. 



(6) 



(7) 



Mfe(o = Ei^^i'' ^ = 1'' 



]<q u\ The g-dcpendence of the so-called generalized fractal di- 

mensions Dq, and therefore a non-linear behaviour of 



TABLE I: The linear system size L, corresponding volume 
V and the number of samples taken and correspondingly the 
total number of wavefunction amplitudes t/j; evaluated. All 
eigenstates are at the critical MIT for which Wc = 16.5 and 
correspond to the five eigenvalues closest to the band center. 
The diagonalization of the matrix was computed on an SGI 
Altix 3700BX2 where for L = 240 it took approximately 24 
hours and requires ~24GB of memory to obtain five eigen- 
states for one disorder realization. 



L V = L samples 



^i 



20 8 X 10^ 24 995 2 x 10** 

30 9 X 10=* 25 025 6.8 x 10*^ 

40 6.4 X lO"* 25 025 1.6 x 10^ 

50 1.3 X 10^ 25 030 3.1 x 10^ 

60 2.2 X 10^ 25 030 5.4 x lO'-* 

70 3.4 X 10^ 24 950 8.6 x 10** 

80 5.1 X 10^ 25 003 1.3 x 10^" 

90 7.3 X 10^ 25 005 1.8 x 10^° 

100 1 X 10" 25 030 2.5 x 10^" 

140 2.7 X 10" 105 2.9 x 10** 

160 4.1 X 10" 125 5.1 X 10** 

180 5.8 X 10" 100 5.8 x 10* 

200 8 X 10" 100 8 X 10** 

210 9.3 X 10" 105 9.7 x 10** 

240 1.4 X lO'^ 95 1.3 x 10^ 



T(g), is an indication of multifractality. Dq is a monoton- 
ically decreasing positive function of q and Dq is equal to 
the dimension of the support of the measure. At critical- 
ity, T{q) can also be parametrized as T(g) = (i((7— 1) + Ag, 
where A^ are the anomalous scaling exponents character- 
izing the critical point.— The singularity spectrum /(a) 
is obtained from the T(g) exponents via a Legendre trans- 
formation, 



dT(g) 
dq ■ 



9 = I' (.a), 



fq = fi^q) = aqq-T{q). 



(8a) 
(8b) 



Here, ,f{a) denotes the fractal dimension of the set of 
points where the wavefunction intensity is IV'iP ^ L^", 
that is in our discrete system the number of such points 
N^ scales as L/("),^ 

The singularity spectrum f{a) is a convex function of 
a and it has its maximum at ao ^ d where /(ao) — 
d. For Qfi we have f{ai) = ai and ,f'{ai) = 1. In 
the limit of vanishing disorder the singularity spectrum 
becomes narrower and eventually converges to one point 
f{d) = d. On the other hand, as the value of disorder 
increases the singularity spectrum broadens and in the 
limit of strong localisation the singularity spectrum tends 
to converge to the points: /(O) = and /(oo) = d. Only 
at the MIT we can have a true multifractal behaviour 
and as a consequence the singularity spectrum must be 
independent of all length scales, such as the system size. 



The symmetry law ([T]) can also be written as 

Qfg -I- ai-q = 2d. (9) 

Since, due to the wave function normalization 
conditiori^^, the singularity strength a is always positive, 
it readily follows that the symmetry requires a ^ 2d. 
Moreover the ^ a ^ d and d ^ a ^ 2d regions of 
the singularity spectrum must be related by f{2d — a) = 
/(a) + d — a, as can be checked by combining Eqs. ^ 
and (191). 



B. Numerical determination of f{a) at the MIT 
using typical average 

The numerical analysis is essentially based on an aver- 
aged form of the scaling law ^ in the limit A = l/L -^ 0. 
This can be achieved either by making the box size I — > 
for a fixed system size L, or by considering L — > oo for a 
fixed box-size. The question of how to compute a proper 
average of the moments Pq is determined by the form of 
their distribution functionj ^^i^^ The scaling law for the 
typical average of the moments Pq is defined as 



,(lnP,(A))^;^r'-(,) 



(10) 



where (■ ■ ■ ) denotes the arithmetic average over all real- 
izations of disorder, i.e. over all different wavefunctions 
at criticality. The scaling exponents are then defined by 

(In A (A)) 



T'^P(q) = lim 



In A 



(11) 



and can be obtained from the slope of the linear fit of 
(InPq) versus In A. Applying Eqs. ([5|) we obtain similar 
definitions for a and f{o:), 



«r = 1™ 1^ /^4('?,A)ln4(l,A) V (12a) 

/fP = lim ^ /x;4(g,A)ln4(g,A)\ , (12b) 

where Skiq, X) = fi1{X)/Pq{X) is the normalized g-th 
power of the integrated probability distribution /Xfc(A). 
The singularity spectrum could also be obtained from 
T{q) by means of the numerical Legendre transformation 
([5]), but this latter method is numerically less stable. 

The typical average is dominated by the behaviour 
of a single (representative) wavefunction. It is because 
of this that the /*^P(a) will usually only have positive 
values, since the average number of points in a single 
wavefunction with a singularity a such that f{a) < 
is L~l-''(")l ^ 1. It is also worth mentioning that due 
to the relation (j8ap . the typical singularity spectrum is 
expected to approach the abscissa axis with an infinite 
slope. However, it has been proven numerically, that the 
region of a values near the ends where the slope tends to 
diverge gets narrower and eventually disappears as the 
thermodynamic limit is approachedi^ 



IV. SCALING WITH BOX SIZE 

In the scaling law of Eq. ([6]), the hniit A ^ can 
be reached by taking the box size I -^ 0, i.e. we are 
evaluating the scaling of Pq with box size I at constant 
L. Numerically, we consider a system with large L and 
we partition it into smaller boxes such that condition 
Im ^ I < L is fuUfilled with Im the lattice spacing. This 
ensures that the multifractal fluctuations of |'0p will be 
properly measured. We usually take values of the box size 
in the range I G [10, L/2]. We have found that the most 
adequate box-partitioning scheme is when the system is 
divided into integer number of cubic boxes, each box with 
linear size ^4^ The system is partitioned in such a way 
that it can be divided equally into boxes and the origin of 
the first box coincides with the origin {x, y, z) — (0, 0, 0) 
of the system. We have used this method to produce all 
the results in this section. We have also tried other box- 
partitioning strategies, however, their results were less 
accurate and will be discussed elsewhere4S 

For each wave function, we compute for the g-th mo- 
ment of the box probability in each box. and Pq as in ([5]), 
as its sum from all the boxes. The scaling behaviour of 
the averaged gIPR with box size pU]) is then obtained by 
varying I. Finally, the corresponding values of the singu- 
larity strength a and spectrum /(a) are derived from the 
linear fits of the Eqs. P^ in terms of the box size. With 
only one system size to be considered, the box-size scaling 
is numerically relatively inexpensive and has been much 
used previously in performing a MFA^ii^iii In Figs. [H 
and [3] we show examples of /(a) and associated linear 
fits. 



A. General features of /*■*''' (a) 

The singularity spectrum for system size L = 240 with 
95 states that is obtained using Eqs. (|12p with I -^ 
is shown in Fig. [2l The /*''? is compared with the cor- 
responding spectrum that is derived from the symmetry 
relation ^ and with the parabolic spectrum^i Here, the 
maximum /'■*'P(ao) = 3 which is equal to the dimension 
of the support can be found very near to ao = 4 where 
the maximum of the parabolic spectrum is located ati^ 
In the region within the vicinity of a = 3, the typical sin- 
gularity spectrum closely resembles the parabolic f{a). 
However, for large \q\ values particularly at the tails, the 
/'^P(a) starts to deviate from the parabolic spectrum. 
We note that the symmetry relation ^ requires that the 
spectrum should be contained below the upper bound of 
a = 2d. 

In order to obtain a and /(a) via the linear fit of 
Eqs. ([T^ . a general x^ minimization is considered tak- 
ing into account the statistical uncertainty of the aver- 
aged right-hand side terms. In this way we can carry 
out a complete analysis of the goodness of the fits via 
the quality-of-fit parameter Q, as well as the usual linear 
correlation coefficient r^ . The behaviour of these quanti- 
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FIG. 2: Singularity spectrum (thick solid black line) obtained 
using box-size scaling of the typical average of Pq for system 
size L = 240 with 95 states. The error bars which are equal to 
one standard deviation mark the locations corresponding to 
integer q values. The corresponding symmetry-transformed 
spectrum ^ is shown as thin black dashed line. The ana- 
lytical parabolic form^ is represented by the thick gray solid 
line. The thin horizontal and vertical lines denote the / = 0, 
/ = 3 and a = 4, a = 6 values, respectively. The values for 
the linear correlation coefficient r^ and quality-of-fit param- 
eter Q for both a'-*''' and /'■*'^(a) are shown in the bottom 
shaded panel. 




FIG. 3: Mass exponents (a) and generalized fractal dimen- 
sions (b) obtained using box-size scaling of the typical aver- 
age of Pq for L — 240 considering 95 states. The filled black 
circles correspond to q = while the empty circles denote 
integer q. The thin horizontal and vertical lines indicate the 
t(0) = —3, Do = 3 and q = values, respectively. In panel 
(c), we show the linear fits of Eq. (|lip for the r*^''(g) in panel 
(a). Only fits for integer values of q ranging from g = — 10 
(top) to q = 10 (bottom) are shown. The values of r'^P(g) 
are given by the slopes of the fits. Data points for g 7^ 
have been properly shifted vertically to ensure optimal visu- 
alization. Data for g = is highlighted with filled symbols. 
Standard deviations of all data are contained within symbol 
size in all cases. 
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FIG. 4; Singularity spectrum obtained using box-size scal- 
ing of the typical average of Pq. In panel (a): system size 
L = 90 for 250 (thick gray line) and 2.5 x 10* (thick black 
lines) states. The corresponding symmetry-transformed spec- 
tra ((21) are shown as thin dashed lines. The error bars mark 
the location of integer q values. In panel (b); system size 
L = 60 (thick gray line), L = 100 (thin gray line), L = 200 
(thin black line) and L = 240 (thick black line) each having 
10^, lO'^, lO'^ and 95 states respectively. The dashed vertical 
line indicates a = 6. In all cases, the error bars denote one 
standard deviation. 



tics for the different parts of the spectrum (correspond- 
ing to different values of the moments q) is shown in the 
bottom panel of Fig. [51 The r^ value is almost equal 
to one for all a which shows the near perfect linear be- 
haviour of the data. Furthermore, acceptable values of 
the Q parameter are also obtained. However, a decline 
in the r^ and Q values is seen at the tails. These re- 
gions correspond to the large \q\ values where the numer- 
ical uncertainties in computing for the Pq over a number 
of different disorder realizations are large enough to af- 
fect the reliability of the data. Figure [31 presents the 
corresponding sets of mass exponents T*yP(g), general- 
ized fractal dimensions D^^p and linear fits for T'^^'^{q) 
for the singularity spectrum in Fig. |21 The g-dependence 
of the decreasing function Dq = T{q)/{q ^ 1) is an indi- 
cation of multifractality. Here, we see that Do — d a.s ex- 
pected. The corresponding T*^^P{q) is shown in Fig.|3l[a). 
It displays the characteristic nonlinearity of a multifrac- 
tal where t(0) = —d. The regions corresponding to large 
\q\ values show a linear behaviour with a constant slope. 
Since the singularity strength is defined as aq = T'{q) 
then a linearity in t'''p((7) found in the limit oi \q\ ^ oo 
results in a'^P values that approach upper a+ and lower 
a- bounds. Hence, the /'^P(a) meets the a^^^ axis at 
these termination points with an infinite slope. Further- 
more, we will show that the location of a_ and a+ is 
greatly affected by system size. For a detailed discussion 
on the relationship between the shapes of r(q) and /(a), 
we refer to the references [ij and l39l. 



B. Effects of the number of states and L 

In panel (a) of Fig. JH we show the small-a region 
of /*yP(a) for the case of L = 90 with 2.5 x 10^ and 
2.5 X 10"' states. When more states are considered for a 
fixed system size, the termination point «_ moves further 
towards smaller values (i.e., towards the thermodynamic 
limit) and the symmetry relation is more closely satisfied. 
However, when a large number of states has already been 
considered (such as 2.5 x 10^ for L = 90) the shape of the 
/*^P(a) will not significantly change anymore with more 
states as illustrated by the already small uncertainties. 
This takes us to consider bigger system sizes in order to 
be able to improve the symmetry relation. In panel (b) 
of Fig.m we show a portion of the large-a part of /^'^^{a) 
for varying system sizes i = 60 with 10"^ states, L = 100 
and L = 200 with 10^ states each, and L = 240 with 
95 states. We see that for the same number of states the 
degree of fluctuations as represented by the size of the er- 
ror bars is larger for smaller system size. Moreover, the 
/'^P(a) for L — 60 with 10^ states is, within the standard 
deviations, the same as that for L ~ 240 with 95 states. 
This can be explained by the total number of wavefunc- 
tion values ■0^ involved in the average, which are nearly 
the same for both cases and hence causes the same shape 
of /*^P(a). This means that when using box-size scaling 
for the typical average of Pq , the number of disorder real- 
izations needed to obtain the singularity spectrum up to 
a given degree of reliability decreases with the size of the 
system. Remarkably, we also see in Fig. IHJb) a general 
tendency that with larger L the singularity spectrum ap- 
proaches the upper bound of a = 6 in keeping with what 
the symmetry relation requires. 

In Fig.|5l we show the spectra corresponding to i = 60 
with 2.5 X 10^ and L ~ 240 with 95 states to clearly 
show the effect of the system size. We observe that the 
value of ao (i-e., location of the maximum) and the shape 
of the singularity spectrum near the maximum do not 
change anymore with L. This L-invariant behaviour of 
the singularity spectrum is an attribute of a critical point. 
In inset (a), with increasing system size the position of 
the termination point a_ moves towards smaller values. 
Furthermore, a closer look of j'^'^^ia) in insets (a) and 
(b) shows that when a bigger system size is used, even 
with less eigenstates, there is a well defined improvement 
to satisfying the symmetry. 



C. Symmetry relation 

In order to quantify how the symmetry is being satis- 
fied with regards to either taking more states or consider- 
ing bigger system size, we present Fig. |6l The top panel 
is an exact calculation of the symmetry relation of Eq. ^ 
whereas the bottom panel shows the difference between 
the singularity spectrum and its symmetry-transformed 




FIG. 5: Singularity spectrum obtained using box-size scal- 
ing of the typical average of Pq, for system sizes L — 60 
(thick gray hue) and L = 240 (thick black hne) with 2.5 x 10'* 
and 95 states respectively. The corresponding symmetry- 
transformed spectra ((2| are shown as thin dashed lines. The 
insets show details for (a) small and (b) large a values. In all 
cases, the error bars are equal to one standard deviation. In 
inset (a), the error bars mark the location corresponding to 
integer q values. 
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FIG. 6; Measures for the degree of symmetry of the multifrac- 
tal spectrum obtained from the typical average with box-size 
scaling. The upper panel shows the numerical evaluation of 
the symmetry law as a function of q for system sizes L = 60 
with 2.5 X 10* states (dashed black), L = 100 with 2.5 x 10* 
(dark grey) and L = 240 (black) with 95 states. The bottom 
panel shows Sf{a) versus a. Note that there is no correspon- 
dence between the abscissa axes of the upper and lower plots. 
For clarity, two values of q for the black line are explicitly 
written. 



counterpart, defined as 



Sfia) = \f{a)~f{2d~a) + d-a\. 



(13) 



The latter plot is an efFective tool to tell us the range 
of the a values where the symmetry is satisfied up to a 
given tolerance. However, directly comparing the degree 
of symmetry via Sf{a) is just an approximation since (i) 
linear interpolation has to be used to measure the verti- 
cal distance properly at several values of a, and (ii) for 
a given q the corresponding value of a as well as its un- 
certainty depend upon sizes and realizations of disorder, 
and this makes the comparison of the different curves in 
terms of a not as reliable as Eq. ^. In fact, the resulting 
error bars are much larger than in the top panel of Fig. [5] 
and even larger than the variation between the 3 shown 
Sf{a) curves. Nevertheless, the results in Fig. [S] illustrate 
that there is a tendency to find a better agreement with 
the symmetry relation whenever more states or bigger 
system sizes are considered. The best situation corre- 
sponds to the biggest system size available (L = 240) 
even though the number of eigenstates is lower than for 
smaller systems. The relatively weak effect of the num- 
ber of states on the shape of the singularity spectrum is 
a result of taking the typical average where by nature 
the average docs not dramatically change with the num- 
ber of samples taken. Furthermore, a rough estimation 
from our results suggest that in order to obtain numeri- 
cally a good /(a) symmetry relation {6f{a) ^ 0.01) for 
a G [1.5,4.5] using box-size scaling one would have to 
consider very big system sizes L ^ 1000. 



V. SCALING WITH SYSTEM SIZE 

The scaling law of the gIPR (O can also be studied 
in terms of the system size L. Obviously the numerical 
calculation of eigenstates for very large 3D systems is a 
demanding taski^i^i^i^ Hence previous MFA studies at 
the MIT have been mostly based on the box-partitioning 
scaling described in Sec. IIVI One naturally would ex- 
pect the scaling with the system size to perform better 
in revealing the properties of the system in the thermo- 
dynamic limit. The fact that for each system size one 
has several independent realizations of the disorder helps 
reduce finite-size effects, which will be unavoidably more 
pronounced when doing scaling with the box size. Ob- 
viously the larger the system sizes and the more realiza- 
tions of the disorder, the better. 



A. Coarse-graining for negative q 

In the present case the scaling variable is L, and the 
formulae (fTTj) and (fT2|) for the singularity spectrum are 
only affected by the substitution: limA^o ^ — linii^oo. 
The box size / which determines the integrated probabil- 
ity distribution fJ.k{l) is now a parameter in the expres- 
sions dm and 1121) for r*yp(g), afP and /fP. Changing 
the value of I is effectively equivalent to renormalize the 
system size to a smaller value L' = L/l. Therefore it 
is clear that the most favourable situation to approach 
the thermodynamic limit is setting 1 = 1, thus defin- 
ing the generalized IPR in terms of the wavefunction it- 



self, Pq = X]i=i iV'iP'- However, when considering neg- 
ative moments, all the possible numerical inaccuracies 
that may exist in the small values of \ipi\'^ will be greatly 
enhanced, which in turn causes a loss of precision in the 
right branch [a > ao) of the singularity spectrum. The 
best way to fix this problem is to use a box-size / > 1 
for g < 0. In this way the relative uncertainties in the 
smallest values of the coarse-grained integrated distribu- 
tion /ifc(/) are reduced with respect to the values of the 
wavefunction. This coarse-graining procedure to evalu- 
ate the negative moments of the wavefunction when do- 
ing system-size scaling was first described in Ref. |23 and 
as we have seen its validity is readily proven when one 
assumes the scaling relation ([6|) as the starting point of 
the MFA. 

The numerical singularity spectrum is thus obtained 
from the slopes of the linear fits in the plots of the av- 
eraged terms in Eqs. (|12p versus InL, for different val- 
ues of the system size L. Where for positive q we have 
/^fe(l) = IV'fcP and for negative q the integrated measure 
/ifc(/ > 1) is kept, with / = 5 in most of the calcula- 
tions. The value of I for the coarse- graining procedure 
should not be very large, otherwise finite-size effects will 
be enhanced again due to the reduction in the effective 
system size. For the benefit of the reader let us rewrite 
the formulae p^ in the particular case where / = 1, 
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for large enough system sizes L. As before the angular 
brackets denote the average over all eigenstates. 



B. General features of /'^^(q) and the effects of the 
number of states and L 

In Fig. [7] we show the singularity spectrum obtained 
from Eqs. HH) {q < 0) and HH) (q ^ 0). We have con- 
sidered system sizes ranging from L = 20 to 100, and 
^ 2.5 X lO'' states for each system size, as shown in Table 
m In spite of the good linear behaviour observed in the 
fits to obtain a'^P and fg^^, shown in Fig.[8l the values for 
Q in the bottom panel of Fig. [71 suggest a loss of reliabil- 
ity near the termination regions of the spectrum. On the 
other hand the standard deviations of the {a'^'P, f^^^{a)} 
values are really small even near the ends. These uncer- 
tainties are directly related to the number of states we 
average over: the more realizations, the smaller these 
uncertainties are. It must be clear that these standard 
deviations only give an idea about the reliability of data 
as a function of the number of disorder realizations for 
the particular range of system sizes that one is using. To 
illustrate the influence of the number of disorder rcaliza- 
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FIG. 7: Singularity spectrum obtained from typical averaging. 
System sizes from L = 20 to 100 have been considered with ~ 
2.5 X 10'' different wavefunctions for each system size, as shown 
in Table ID The dashed line is the symmetry-transformed 
spectrum according to /(6 — a) = f{a) + 3 — a. The values 
of q range from g = — 10 to q = 10 with a step of 0.1 {1 = 1 
for q ^ and I = 5 for q < 0). Error bars which are equal 
to one standard deviation highlight the values corresponding 
to integer q. The lower panel shows the linear correlation 
coefficient (r^) and the quality-of-fit parameter (Q) of the 
linear fits to obtain the values for a and f{a). 



tions upon the typical average a comparison can be found 
in Fig. [HI between the P^^{a) spectrum obtained after 
averaging over 10'^ states for each system size and the one 
for ^ 2.5 X lO'' states. As can be seen, after this increase 
in the number of states the overall change in the spec- 
trum is not very significant, altough some variation can 
be noticed in the regions shown. In particular, the right 
branch of the spectrum moves inwards and the end of 
the left tail shifts to smaller values of a. In both regions 
the expected variation of the spectrum is well described 
by the standard deviations. In the case of Fig. [7] accord- 
ing to the standard deviations the conclusion is that a 
further increase of the number of states will not mean a 
significant change in the shape of the spectrum. Never- 
theless it must also be very clear that if we consider a 
different range of larger system sizes, noticeable changes 
could happen in the singularity spectrum. The standard 
deviations do never account for the effects stemming from 
the range of system sizes used. 

To evaluate the effects due to the system size, we com- 
pare in Fig. [TO] the multifractal spectrum obtained con- 
sidering different ranges of system sizes with a similar 
number of disordered realizations. In the main plot it 
can be seen how the shape of the spectrum changes in its 
right (large a) branch, which moves inwards, when we 
consider system sizes in the interval [140, 240] compared 
to the situation for sizes in [40, 100]. The left end of the 
spectrum also shifts to smaller values of a when larger 
system sizes are considered. In this case the standard 
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FIG. 8: Linear fits of Eqs. ^ for af p values (left) and /f p 
values (right) of the singularity spectrum in Fig. [T] Only fits 
for integers values of q ranging from q = 10 (top) to g = — 10 
(bottom) are shown. The values of a^^P and /^-^'p are given by 
the slopes of the fits. Data points for g 7^ have been properly 
shifted vertically to ensure optimal visualization. Data for 
q = highlighted with filled symbols. Standard deviations 
are contained within symbol size in all cases. 
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FIG. 9: Left branch (a) and right branch (b) of the singularity 
spectrum obtained from typical average, scaling with system 
sizes from L = 20 to 100 for ~ 2.5 x 10"* states (black) and 
10^ (grey) for each system size. Dashed lines correspond to 
spectra transformed according to the symmetry law. The 
values of q range from g = — 10 to q = 10 with a step of 0.1 
(/ = 1 for q ^ and I = 5 for q < 0). Error bars are standard 
deviations. 



deviations are noticeable since we have only considered 
10^ states for each system size. In the insets (c) and 
(d) within Fig. [10] a similar comparison can be found for 
ranges of smaller sizes, [60, 100] versus [20, 60] but with 
a much higher number of states, ^ 2.5 x lO** for each 
size. In this situation the change is less dramatic, but 
the tendency remains the same. In particular it should 
be noticed in Fig. [TUf c) how the change in the left end 
of the spectrum is not contained in the uncertainty re- 
gions given by the error bars, confirming the fact that 




FIG. 10: Singularity spectrum obtained from typical average 
using different ranges of system sizes. Grey line: 7 system 
sizes from L = 40 to 100 and 10 states for each. Black line: 
6 system sizes from L = 140 to 240 and ~ 10^ states for each. 
Insets (c) and (d): Grey line: 5 system sizes from i = 20 to 
60 and ~ 2.5 x 10* states for each. Black line: 5 system sizes 
from L = 60 to 100 and ~ 2.5 x 10** states for each. In all 
cases dashed lines correspond to spectra transformed via the 
symmetry law. The values of q range from g = — 10 to 5 = 10 
with a step of 0.1 {1 = 1 for g ^ and I = 5 for g < 0). Error 
bars are standard deviations. 



these standard deviations do not fully describe system 
size effects. 



C. Symmetry relation 

The symmetry relation ^ is only partially fulfilled in 
Fig. [71 Still, a nice overlap between the original spectrum 
and the symmetry-transformed one occurs in the region 
around the symmetry point a = 3. The agreement is lost 
when approaching the tails, which are the parts more af- 
fected by numerical inaccuracies and system-size effects. 
For a given range of system sizes, the symmetry relation 
tends to be better satisfied whenever the number of dis- 
ordered realizations is increased, as can be observed in 
Fig. [21 On the other hand the improvement of the sym- 
metry is even more dramatic when we consider larger 
system sizes to do the scaling, as shown in the insets 
(a) and (b) of Fig. [TUj In this figure it is evident how 
the value of /(a = 6) decreases when considering larger 
system sizes, hence tending towards the upper bound at 
a = 2d as predicted by ([2]). 

A quantitative analysis of the symmetry relation is 
shown in the upper panel of Fig. 1111 The best data corre- 
spond to the scaling with system sizes in [140, 240] after 
averaging over ~ 100 states for each size (cp. Table [I|). 
Even with such a low number of disorder realizations, the 
observed symmetry is better on average than the one ob- 
tained for sizes in [20, 100] with 2.5 x 10^ states for each 
L. Let us emphasize that for L — 100 the total num- 
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FIG. 12: Comparison of degree of symmetry for tfie spectra 
obtained using typical averaging of Pq for the cases of box- 
size [BS] scaling (gray) and system-size [SS] scaling (black). 
The best spectrum for each case has been considered: [BS] 
L = 240 with 95 states and [SS] L G [140,240] with 100 
states for each size. The plot shows the numerical evaluation 
of the symmetry law as a function of q. 



FIG. 11: Measure of degree of symmetry of the multifractal 
spectrum of Fig. [10] obtained from scaling with system size. 
The upper panel shows the numerical evaluation of the sym- 
metry law as a function of q. The bottom panel contains 
Sf{a) versus a. Dashed black: 7 system sizes from L = 40 to 
100 and 10^ states for each. Grey: 9 system sizes from L = 20 
to 100 and ~ 2.5 x 10* states for each. Solid black: 6 system 
size from L — 140 to 240 and ~ 10^ states for each. There 
is no correspondence between the abscissa axes of the upper 
and lower plots. For clarity, two values of q for the black line 
are explicitly written. 



ber of wave function values involved in the calculation is 
2.5 X 10^° while for L = 240 it is only 1.3 x 10^. This 
shows that although the number of disorder realizations 
is important to improve the reliability of data (reducing 
the standard deviations), the effect of the range of system 
sizes is more significant. And although it can be argued 
that the error bars of the black line in the upper panel 
of Fig. [11] are still very large, we have already shown 
that when increasing the number of states the symme- 
try simply gets better and thus the line will move even 
closer to zero. In the lower panel of Fig. [11] the deviation 
from symmetry Sf{a) defined in (|13p is also shown and 
corroborates these findings. 

Hence, whenever the reliability of data is improved by 
increasing the number of disorder realizations, or when 
finite-size effects are reduced by considering larger system 
sizes, we get a better agreement with the symmetry law 
([2]) of the multifractal spectrum. Assuming the degree of 
symmetry is a qualitative measure of the MFA itself, then 
from a numerical viewpoint, the best strategy when doing 
scaling with system size and typical averaging would be 
to go for the largest system sizes accessible even though 
it means having less realizations of disorder. 



disorder realizations or by taking larger system size, the 
/*^p(q;) spectrum becomes evermore close to obeying the 
proposed symmetry relation ([2]). Using the typical av- 
erage, the best symmetry in the singularity spectrum is 
obtained by taking large system sizes. Due to the nature 
of the typical averaging, taking more states only changes 
the shape of the /'''^(a) up to a point. By consider- 
ing larger system sizes, a significant improvement of the 
symmetry relation is achieved, leading to lower values of 
q;*^p and /'^p on the left side of the spectrum as well as 
a better agreement with the upper cut-off of a < 6. 

In Fig. [121 let us now compare box- and system- 
size scaling. With system-size scaling the symmetry is 
(nearly) satisfied for a wider range of a values as com- 
pared with the box-size scaling. Box-size scaling is more 
strongly influenced by finite-size effects. However, the 
agreement with the symmetry relation is lost for both 
methods at large jgj or equivalently at |a — 3| ^ 0. Un- 
surprisingly, these are the regions greatly affected by nu- 
merical inaccuracies and finite-size effects. Hence we con- 
clude that within the accuracy of the present calculation 
and within the limits of the typical averaging procedure, 
the proposed symmetry relation ^ is valid at the An- 
derson transition in 3D.— 

Last, let us remark that the relation ^ implies neg- 
ative values of / for small values of a. As discussed 
previously, this is hard to see using the typical averaging 
procedure. In Ref. [2^, we have also performed MFA us- 
ing the ensem6/e-averaged box- and system-size scaling 
approaches. The results again support the existence of 
the symmetry ([ij for an even larger range of a values 
and including a negative /(a) part for small a. 
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